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We use the recently introduced small-world networks (SWN) to 
model cross-linked polymers, as an extension of the linear Rouse- chain. 
We study the SWN- dynamics under the influence of external forces. 
Our focus is on the structurally and thermally averaged SWN stretch- 
ing, which we determine both numerically and analytically using a 
psudo-gap ansatz for the SWN-density of states. The SWN stretch- 
ing is related to the probability of a random-walker to return to its 
origin on the SWN. We compare our results to the corresponding ones 
for Cayley trees. 
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I. INTRODUCTION 



Topological properties of polymers can dramatically effect their dynamical properties, 
such as their collapse in bad solvents and their response to external force^Hi. Such 
forces can be applied microscopically, either by having charged polymers (polyelectrolytes, 
polyampholytes) in electrical fields, or via optical tweezers or magnetic bead^l^. In 
this communication we study the stretching of cross-linked objects, whose backbones are 
regular lattices (we will consider for simplicity a ring), a few elements of the backbone 
being chemically connected to each other via cross-links. An experimental realization may 
be a very dilute solution of linear chains which are then cross-linked by irradiation! or 
chemically. Our model is a realization of the so-called small-world networks0Hii (SWN), 
and the disorder (the cross-links) is in statistical terms quenched. Considering all bonds 
to be equal we study the dynamics in the framework of the Rouse modelffl, in which 
the monomers are connected by harmonic springs; we term this structure the small- 
world Rouse network (SWRN). Our SWRN is built out of an N-monomer ring with 
superimposed fixed links between randomly chosen monomers; no additional links are 
generated or broken in external fields. The SWRN is a new intermediate between linear 
chains and networks. Distinct from Cayley-trees, which model hyperbranched polymers 
without loops, loops are a fundamental ingredient here. As such the SWRN is interesting 
in its own right as a study of the interplay between dynamics and topology, and it belongs 
to the class of generalized Gaussian structuresEHi. 



2 



II. DYNAMICS OF SMALL- WORLD ROUSE NETWORKS 

The construction of the SWN which we consider here (see also§), is shghtly different 
from the original one by Watts and Strogatz0, but preserves its main SWN-characteristics. 
Starting from a ring of N monomers, we cross-hnk with probabihty p each monomer 
randomly to any of the monomers of the network. Such cross-links thus connect monomers 
far apart along the chemical backbone (the ring), rendering them close in Eucledian 
space. The SWRN-monomers are, in accordance with the Rouse modelffl, connected 
by harmonic springs of strength k. The position Unit) of the nth bead under the action 
of an external force Fn{t) and in the presence of thermal noise T]^{t) is governed by the 
Langevin equation 

= fc5:A„,R,(t) + F„(t) + r7„(t). (1) 

Here 7 is the coefficient of friction, and the matrix Aij is the connectivity matrix of the 
network. It is defined as follows: Every connection between site i and site j contributes 
—1 to Aij, while An and Ajj are determined from the condition that J2j ^ij = J2i ^ij = 
0. Monassonlii recently published a detailed study of the spectrum of the small-world 
connectivity-matrix A. Among his findings was the existence of a "pseudo-gap" in the 
SWN density of states p{E) which has the form 

C \ 



p{E) ~ 



E^O (2) 



EJ 

This behavior could not be confirmed numerically through direct diagonlization, and 
appeared in the data as a real gap. There are, as we shall see, other ways of probing this 
behavior numerically, and our results support that p{E) behaves as Eq. (0). 
Writing 



R(t) = (Ri(t),R2(t),...,R^(t)f (3) 

and similary for the other quantities in Eq. (|1]), we can rewrite the equation of motion in 
a more condensed form as 

dK(t) F(t) 
at 7 

Here we have introduced a = k/j and w = 77/7. In the case of a spatially constant 
external force, the solution to Eq. (H) is obtained as 

R(t)=y'* dse-''^^'-'^(^ + w{s)^ (5) 

We now specialize to the following situation: the force is switched on at time t = and 
pulls only the mth bead in the ?/-direction, i.e. Fj(t) = 9(t)Si^rnF y. Here 9(t) is the 
Heaviside step-function, 6ij Kronecker's delta and y is a unit-vector pointing in the y- 
direction. We focus on the displacement of the mth bead along the y-axis, and average 
over thermal noise, using (w(t)) = 0. Finally we perform a structural average over m and 
end up with (see e.g.i^ and references therin for details) 

In this equation Rm,y is the y component of R^ and Aj with i = 1 . . . N are the eigenvalues 
of the connectivity matrix A. The last expression follows from the fact that for a connected 
structure, only one eigenvalue vanishes (say, Ai). At times t much smaller than the time 
scale set by the largest eigenvalue Xmax, i-e. when aXmaxt ^ 1, y(t) increases linearly in 
time: Y{t) ~ Ft/'j. That is, only the monomer being pulled moves with a constant speed, 
not yet feeling the influence of the other monomers. Likewise, at late times aXmmt ^ 1, 
where Xmin is the lowest non-vanishing eigenvalue, the entire polymer is being pulled with 
a constant speed, Y{t) ~ Ft/{N'j). These observations are independent of the specific 

4 



structure being pulled, and only in the intermediary regime X^ax ^ crt <^ A~^„ does the 
particular topology of the polymer affect the dynamics, namely through the spectrum of 
the connectivity matrix. 

The numerical computation of the quantity Y{t) above proceeds as follows. From a 
specific realization of an = 1000 small-world network, we construct the correspond- 
ing connectivity matrix. Then we find the N eigenvalues using standard routines, and 
implement Eq. (^). To get an idea of the importance of sample to sample fluctuations, 
we consider first 10 different realisations of the SWRN for p = 0.05. Plotted is in Fig. 
m on double logarithmic scales y as a function of t, where here and in the following we 
use the dimensionless variables Y*{t) = a'~fY{t)/F and t* = at. In Fig. |I| we display the 
envelope of all 10 realizations, i.e. the two curves are the extremal two "worst" cases. We 
see that the difference in the results is quite small (and appears, as it should, only for 
intermediary times), and therefore we regard results from any specific realization as being 
typical. In Fig. ^ we analyze the dependence of Y{t) on cross-linking, by varying p from 
p = (standard Rouse-model of the ring) to p = 0.01 and p = 0.05. We note first that the 
differences are now considerably larger than in Fig. |I|. Second, for p = i.e. for the Rouse 
chain, we have the standard picture: a subdiffusive \/t behavior at intermediary times 
is followed by a diffusive t behavior at longer times. At very early times we also have a 
linear behavior, albeit not visible in the range of the figure. The initial and final dynamics 
are in accordance with the explanation given abovei. The intermediate behavior reflects 
the structure of the spectrum of A, and is also well understood in the Rouse case; it is a 
result of the rather slow propagation of disturbances through the chain (here the ring). 

We infer from the other curves in Fig. |^, that even very small but nonvanishing p 
affect the intermediate behavior of Y{t) quite strongly. For increasing p the curves bend 
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downwards from the p = case, mirroring the increased stiffness of the polymer due to the 
additional links. Thus a ring with cross-links can be easily distinguished experimentally 
(say through NMR or electronic energy transfer) from one without cross-links, whose Y{t) 
dynamics under F is never slower than ^/t. As expected, the very early and very late 
behavior in all three cases coincide, being independent of the specific structure under 
scrutiny. 

In Fig. ^ we plot on logarithmic scales Y*{t) as a function of time t* for several values 
of p. Increasing p increases the stiffness of the polymer, and this is reflected in the 
intermediary regime, which becomes almost fiat for large p. Moreover the long time 
behavior of the polymer is reached much earlier for polymers with a higher number of 
cross-links. In line with the discussion above of the range of the intermediary regime 
i^max ^ crt <^ -^mL)) ^his feature means that the lowest non- vanishing eigenvalue Xmin 
gets quite large, and it is hence related to the appearance of a (pseudo) gap in the spectrum 
of A. 

III. THEORETICAL ANALYSIS 

As indicated earlier, the initial as well as the asymptotic behavior of Y(t) are well 
understood. Thus we will concentrate here on the richer and much more complex in- 
termediate behavior. We shall rewrite Eq. (j^) using a continuous picture, based on the 
density of states p(A) = lim^^oo{^/N) — K), but continue to separate out the 

vanishing eigenvalue from the rest. Hence, with e very small, e O"*": 

F Ft F /■* 

ym-J^l''.Ze-''"^- + -ldsl <;Ap(A)e-- (7) 

an expression which is a fortiori correct in the presence of a gap, where one can take 
< e < Xmin- It will be convenient also to consider the stretching (relative motion) A{t) 
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separately: 

Ft F /■* 
A{t)^Yit)- — = - ds d\p{X)e 

7 JO Jo+ 

We remark that the inner integral in Eq. @ is related to the probability for a random 
walker to be present at the original siteE'El. We therefore first analyse the behavior of 
this quantity: 

POO 

Po{t)= / dXp{X)e-^' (9) 



The asymptotic temporal behavior is accessed through the behavior of p(A) at small A. 
Inserting the expression of Monassonlll Eq. (0) into Eq. (^) we obtain: 
P„W^/;«A-V^exp(-^-«) ^-r=/=A^~,,,.pp3(^_^)) (10) 

The asymptotic behavior of the integral follows readily from a saddle-point procedurelll, 
so that we end up with 

Po{t)r^t-'/'expi-C't'/'), (11) 

whith C = 3 (C/2)^^^.This is by itself a quite interesting and novel result, and it compares 
favourably to our numerical simulationsO for Po{t). The result is close in form to that 
for Cayley trees0il, where Po{t) ~ exp(— ct). Inserting Eq. in Eq. (||) and 

reintroducing a, we get 

Notice that for very large t the stretching A{t) of the SWRN tends to a constant Aqo- In 
the units of our figure this constant depends mainly on C (since F, 7 and a drop out); 
theoretically one may obtain C and also C out of Aqo- 

In Fig. ^ we plot the dimensionless stretching A*(t) = Y*{t) — at/N for p = 0.05 and 
compare it to the analytical form Eq. ([12|). We do this by fitting a — bt^^^ exp(—ct^^^) 



to the data, and as can be inferred from Fig. ^, the agreement is very convincing. From 
the least-squares fit we obtain a = 5.09, b = 7.67 and c = 0.54. We remark that the 
agreement at short times may be rendered even better by also keeping the next term in 
the expansion of the integral of Po{t), a term which is proportional to t^^^^ exp (^—C't^^^^ . 
Furthermore we remark that for Cayley-treesS the intermediate behavior of Y(t) can also 
be determined in a similar manner: The saddle-point approach yields to leading order: 

Y{t) ~ a - br^/^ exp(-at), (13) 

with a, b and c being constants. 

IV. CONCLUSIONS 

In this communication we have studied the behavior of a small-world network model 
(the SWRN) of a linked ring-polymer, focusing on its dynamics under external forces. 
Thus the motion of a monomer pulled by such a force is vastly different, depending on 
whether the monomer belongs to a SWRN or to a simple ring without cross links. This 
may enable via NMR or electronic energy transferll to distinguish clearly between cross- 
linked and non-cross-linked polymers. Our numerical results for the stretching of the 
SWRN under external forces are in excellent agreement with our analytical expressions, 
which used the pseudo-gap behavior of the SWN density of states, as postulated in former 
wor k§. As discussed in the present communication, these results are directly connected 
to expressions for the return to the origin of a random walker on the SWN. 
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CAPTIONS 



FIG. 1. Two different realizations give rise to similar behavior of Y{t), here plotted on 
logarithmic-logarithmic scales for p = 0.05. 

FIG. 2. On double logarithmic scales we plot the position Y{t) as a function of time. 
From upper to lower curve, p — 0, p — 0.01 and p — 0.05. 

FIG. 3. On double logarithmic scales we plot the position Y{t) as a function of time. 
From upper to lower curve, p — 0.01, p — 0.05, p — 0.1, p — 0.2, p — 0.5 and p — 0.8. 

FIG. 4. Comparison of the theoretical prediction (dashed) with the data (dash-dotted), 
for p = 0.05. 
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